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ABSTRACT 

Hydro-dynamical modeling of the inner Galaxy suggest that the radius of the outer Lind- 
blad resonance (OLR) of the Galactic bar lies in the vicinity of the Sun. How does this 
resonance affect the distribution function in the outer parts of a barred disk, and can 
we identify any effect of the resonance in the velocity distribution actually observed in 
the solar neighborhood? To answer these questions, detailed simulations of the velocity 
distribution, f(v), in the outer parts of an exponential stellar disk with nearly flat rota- 
tion curve and a rotating central bar have been performed. For a model resembling the 
old stellar disk, the OLR causes a distinct feature in f(v) over a significant fraction of 
the outer disk. For positions up to 2kpc outside the OLR radius and at bar angles of 
~10-70 degrees, this feature takes the form of a bi-modality between the dominant mode 
of low- velocity stars centred on the local standard of rest (LSR) and a secondary mode of 
stars predominantly moving outward and rotating more slowly than the LSR. 

Such a bi-modality is indeed present in f(v) inferred from the Hipparcos data for 
late-type stars in the solar neighborhood. If one interpretes this observed bi-modality as 
induced by the OLR - and there are hardly any viable alternatives - then one is forced 
to deduce that the OLR radius is slightly smaller than R . Moreover, by a quantitative 
comparison of the observed with the simulated distributions one finds that the pattern 
speed of the bar is 1.85 ± 0.15 times the local circular frequency, where the error is 
dominated by the uncertainty in bar angle and local circular speed. 

Also other, less prominent but still significant, features in the observed f(v) resemble 
properties of the simulated velocity distributions, in particular a ripple caused by orbits 
trapped in the outer 1:1 resonance. 

Subject headings: Galaxy: kinematics and dynamics - - Galaxy: structure - - solar 
neighborhod 
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1. Introduction 

The fact that the Milky Way is barred has been estab- 
lished first of all by the interpretation of the gas velocities 



observed in the inner Galaxy ( ie Vaucouleurs 196 z 




ters 1975; Cohen & Few 1976|; Liszt & Burton 198C ; 


Gcr- 


lard & Vietri 1986|; Mulder & Liem 1986; Binney 


et al. 


1991) and later confirmed by infrared 
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k Spergel 1991; Weiland et al. 1994; 


Dwck et al. 1995; 


Binney Gerhard & Spergel 1997) and asymmetries in the 



distribution or magnitude of stars ( 
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Weinberg 1992 




Nikolaev 


Sz Weinberg 1997; 


Stanek 1995 


Sevenster 1996 


; 3tanek 


et al. 1997 


). However, there is still substantial debate on 



the structure and morphology of the bar, its orientation 
with respect to the Sun, and the rotation rate, usually 
referred to as pattern speed. Recent hydro-dynamical in- 
vestigations of Englmaier & Gerhard (1999) and Weiner 
& Sellwood (1999) and the combined stellar- and gas- 
dynamical models of Fux (1999a) suggest that the bar 
rotates fast, i.e. that corotation occurs somewhere be- 
tween 3.5 and 5kpc (for Rq = 8kpc) not much beyond 
the end of the bar. Moreover, the bar angle, the azimuth 
of the Sun with respect to the bar's major axis, is re- 
stricted both by photometry and kinematics to lie in the 
range between about 10° and 45°. 

Apart from the corotation resonance (CR), a rotating 
bi-symmetric bar creates two other fundamental reso- 
nances, the inner (ILR) and outer (OLR) Lindblad res- 
onances, which occur when 



(1) 



Here, u)^ and uir are the azimuthal and radial orbital 
frequencies, while fib is the pattern speed of the bar. An 
orbit for which relation ([l]) is satisfied is in phase with 
the bar, i.e. after each completion of one radial lobe the 
orbit is at the same position with respect to the bar. 
Thus, a star on a resonant orbit is always pushed in the 
same direction, and thus forced into another orbit. 

For circular orbits, cj^ = fl(R) (circular frequency) 
and ujr — k(R) (epicycle frequency) and each resonance 
can be mapped uniquely to the radius of the resonant 
circular orbit; hereafter R lr and R CR denote the radii 
where circular orbits are in outer Lindblad and corota- 
tion resonance. For a nearly flat Galactic rotation curve, 
]J OLR w 1.7i? CR and, from the aforementioned estimates 
for R C n, we find -R lr ~ 6-9 kpc. Thus, we expect the 
OLR of the Galactic bar to be located in our immediate 
Galactic surrounding. The goal of this paper is to an- 
swer the question whether and how this proximity of the 
resonance affects the stellar kinematics observable in the 
solar neighborhood. 




Fig. 1. — Closed orbits (solid) just inside and outside the OLR 
of a rotating central bar (shaded ellipse). The circles (dotted) 
depict the positions of the ILR, CR, and OLR (from inside out) 
for circular orbits. Note the change of the orbits' orientation at the 
OLR, resulting in the crossing of closed orbits at four azimuths. A 
possible position of the Sun is shown as filled circle. The bar angle 
<j> is indicated for the case of a clockwise rotating bar. 

Stellar Dynamics Near the OLR 

Using linear perturbation theory for near-circular orbits 
(cf. Binney & Tremaine 1987, p. 146-151), one finds that 
closed orbits in a barred potential are elongated either 
parallel or perpendicular to the bar. The orientation 
changes at each of the fundamental resonances: inside 
ILR orbits are anti-aligned (so-called X2 orbits), between 
the ILR and CR they are aligned (xi orbits), while be- 
tween CR and OLR the orbits are anti-aligned, until they 
align again beyond the OLR. 

The situation at the OLR is sketched in Figure [l], 
which shows two closed orbits (solid curves) just in- and 
outside the OLR as they appear in a frame of refer- 
ence corotating with the bar. In this frame, the orbits 
near OLR rotate counter-clockwise for a clockwise ro- 
tating bar, like that of the Milky Way. Thus, at bar 
angles <f> between 0° and 90°, the closed orbits inside 
OLR move slightly outwards, while those outside OLR 
move inwards. Clearly, if all disk stars moved on closed 
orbits, the stellar kinematics would deviate from that of 
a non-barred galaxy only at positions very close to R i.m 
where the closed orbits are significantly non-circular. In 
particular, at azimuths where the closed orbits from ei- 
ther side of the OLR cross, one would expect two stellar 
streams, one moving inwards and the other outwards^]. 

However, in general stellar orbits are not closed, but 
exhibit radial oscillations. Many of these orbits are 



trapped into resonance (Weinberg 1994), and may be de- 



1 Based on this consideration, Kalnajs (1991) suggested that the 
Hyades and Sirius stellar streams in the solar neighborhood are 
caused by the Sun being located almost exactly at one of the two 
possible positions in the Galaxy where such orbit crossing occurs 
with the appropriate sign of the radial velocity difference. 
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scribed, to lowest order in their eccentricity, as epicyclic 
oscillations around a closed parent orbit. This means in 
particular, that such trapped eccentric orbits from inside 
the OLR, i.e. with fit, < ui^ + ^wr, can visit locations 
outside i? LRi and vice versa. Thus, the stellar velocity 
distribution observable in the solar neighborhood, fo(v), 
shall be affected by the OLR whenever near-resonant or- 
bits pass near the Sun and are sufficiently populated. 
We may estimate, using the epicycle approximation in a 
back-of-the-envelope calculation, that for late-type disk 
stars this condition is satisfied for -R lr in the range 
fromf] 6 to 9kpc (for i?n = 8kpc). This coincides with 
the above estimate for i? LR- Thus, we expect the OLR 
of the Galactic bar to affect the velocity distribution ob- 
servable locally for late-type stars. Based on the proper- 
ties of the closed orbits, one would expect also for stars 
moving on non-closed orbits different typical velocities 
depending whether they originate from in- or outside the 
OLR. Thus, the expected effect on fo(v) is a bi-modality. 

Weinberg (1994) has studied the orbital response to 
a rotating bar using epicycle theory and assuming a flat 
rotation curve. He found indeed that different orbital 
families may overlap creating a bi-modality in the veloc- 
ity distribution (though he only considered the resulting 
increase in the velocity dispersions) . Weinberg also con- 
sidered a slowly decreasing pattern speed, which affects 
the relative number of stars trapped into resonance with 
the bar, but leaves the final phase-space position of the 
resonance, and hence the velocity of possible modes in 
fo{v), unchanged. 

In order to verify and quantify these estimates and 
expectations of the response of a warm stellar disk to a 
stirring bar and its OLR, I performed numerical simula- 
tions, presented in Section g. In Section ||, I discuss the 
closed orbits in the outer part of a barred disk galaxy 
and their relation to the features apparent in the local 
f(v), while Section ^ gives quantative estimates for the 
velocity of the secondary mode induced by the OLR. In 
Section ||, the velocity distribution observed in the so- 
lar neighborhood, which indeed shows a bi-modality, is 
compared both qualitatively and quantitatively to those 
emerging from the simulations. Finally, Section ^ con- 
cludes and sums up. 

Throughout this paper, I will use units of kpc and 



2 The RMS epicycle amplitude of stars in a population with radial 
velocity dispersion <r fi is X pa v / 2o".fl/ K - Stars originating from a 
radius R = Rq + x could visit us, if \x\ < X(R) yielding 



\x\ < V2 



k(Rq) 



x R~ 



■n- 



where I have assumed that crji decays exponentially with scale 
length Ra and that k <x R _1 (as for flat rotation curves). With 
o-r(Ro) SSkms" 1 , ft(iJo) ~ 35 km s _1 kpc" 1 , and -Ro = Rrr = 
8 kpc, this results in the two solutions x — 2 and 1 kpc. 



kms 1 for radii and velocities, while frequencies and 
proper motions are given in kms -1 kpc -1 . 

2. The Simulations 

In order to investigate the typical structure of a dynam- 
ically warm stellar disk in the presence of a rotating 
bar, we simulate the slow growth of a central bar with 
constant rotation frequency in an initially axisymmetric 
equilibrium representing a warm stellar disk. Clearly, 
this does not simulate the true evolution of the Galaxy, 
since the stellar disk was hardly dynamically warm al- 
ready when the bar formed. Moreover, the bar has pre- 
sumably developed too, both in strength and in pattern 
speed. However, we are not aiming at simulating the 
formation of the Milky Way, but at answering the ques- 
tion for the effect of the OLR, which significantly alters 
the phase-space structure, independently of the forma- 
tion history of the Galaxy. 

2.1. The Simulation Technique 

Currently, traditional iV-body simulations are not suit- 
able for studying the influence of the central bar on a 
stellar disk at high resolution: the required number of 
particles exceeds 10 8 already for merely two-dimensional 
simulations, which makes them very CPU-time intensive 
and aggravates any investigation of the parameter space. 
However, iV-body simulations are unnecessary luxury, as 
they provide dynamical self-consistency and model the 
whole system, both of which is not required in this study. 
Instead, we (i) assume the growth of a central bar, rather 
than simulate a bar-instability self-consistently, and (ii) 
only integrate trajectories which are crucial to the lo- 
cal velocity distribution, rather than those for the whole 
system. Moreover, Poisson noise is avoided by placing 
the trajectories on a regular grid in the observables and 
integrate backward in time (see Section 2.1.3 for details). 

This technique, which may be called backward in- 
tegrating restricted N-body method, is similar in spirit 
to that used in the method of perturbation particles 
( Leeuwin, Combes fc Binncy 1993|) , and has first been 



used by Vauterin & Dejonghe (1997, 1998), who studied 
the non-linear evolution of the stellar distribution func- 
tion (DF) in the inner parts of a bar-unstable model, 
whereby using the analytic potential due to the domi- 
nant linear bar mode. 

Because we are only interested in the effects on the 
planar motions in the outer disk, our modeling is two- 
dimensional, and we do not care much about the details 
of the inner Galaxy. 

2.1.1. The Initial Equilibrium 

Since we are interested only in near-circular trajectories 
passing through a point in the outer parts of the disk, 
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a possible inadequateness of our model potential in de- 
scribing the inner parts of the Galaxy is unimportant. 
Therefore, we take a simple power-law rotation curve 



v c = v (R/Rof 



(2a) 



as our initial axisymmetric model. Here, Ro denotes the 
Sun's distance from the Galactic center and vo the local 
circular speed. The corresponding underlying potential, 



MR) 



Vn X 



f (2/3)- 1 (R/R ) 2f3 
\ ln(R/R ) 



for (3 ± 
for (3 = 0, 



(2b) 



is meant to originate partly from the gaseous and stellar 
components (disk, bulge/bar, halo) and partly from dark 
matter. In restricting our model to two dimensions, we 
make no assumptions about the vertical distribution of 
the matter generating the potential (2b). 

The dynamical equilibrium of the stellar disk is de- 
scribed by a simple DF, / = F(E,L Z ), which depends 
on energy E and angular momentum L z of the stars, 
and is given by equation (10) of Dehnen (1999b). The 
parameters of the DF are set such that it generates, to 
very good approximation, an exponential disk with scale 
length R s and an exponential velocity-dispersion profile 
with scale length R a and local radial velocity dispersion 
0o = &r(Ro)- 

2.1.2. The Bar Potential 

The contribution of the bar to the potential in the outer 
disk is dominated by the quadrupole^], and we neglect 
higher-order multipoles for simplicity. Some orbits vis- 
iting the outer disk may have passed through the inner 
disk and bar region. Therefore, we will use a slightly 
more elaborate model for the bar potential than that 
resulting from just a constant quadrupole moment: 

[-{R h /Rf HR>R h , 
4> h = A h cos(2[4 ) -n h t})x{ (3) 

{(R/R h ) 3 -2 iiR<R h , 

where Rb and Sib are the size and pattern speed of the 
bar. The amplitude A\> of the quadrupole is switched on 
smoothly. It is zero before t — 0, grows with time at 
< t < t\ as 

A h =A i (—e-~e+—ti+-), c = 2— - 1, (4) 

and stays constant at ^4b = Af after t\. Thus, the am- 
plitude and its first and second time derivative behave 
continuously for all t, guaranteeing a smooth transition 
from the non-barred to the barred state. 



3 Bar formation is mainly a re-arrangement of the matter inside 
corotation such that the monopole of the potential outside this 
region remains unchanged. 



2.1.3. The Time Integration 

The collisionless Boltzmann equation tells us that the 
DF remains constant along stellar trajectories. Thus, 
the value f(w,t2) of the DF at some phase-space point 
w = (x, v) and time t% is equal to f{w , 0) if w originates 
from integrating Wq from t = to t = t%. In our case 
and in contrast to iV-body simulations, the potential as 
function of space and time is known a priori. Therefore, 
we can obtain Wq just by integrating the orbit passing 
through w at ti backward in time until t = 0. 

In practice, we specify the final time ti and choose 
the phase-space points w from a grid in planar velocities 
at a position (Ro,<fi) in the Galactic plane. Each of the 
resulting orbits is integrated backward until t = 0, and 
the initial energy E and angular momentum L z are re- 
membered. From these, the value /(«?, t%) = f(wo,0) = 
F(E, L z ) can be computed for any initially axisymmetric 
equilibrium DF F. 

2.2. The Parameters 

Some of the parameters arising in the simulations are 
unimportant for our purposes. For instance, the size 
i?b of the bar, which, in agreement with the findings 
mentioned in Section [l], we fix to be 80% of the corotation 
radius. For our model, the latter is given by 



i? CR EEi? (r» b /r!o) 1/(/3 " 1) , 



(5) 



where fio =vo/ Ro denotes the local circular frequency. 

Since our simulations are not self-consistent, the pa- 
rameters of the DF are unimportant for the dynamics of 
the stellar orbits, but crucial for their population with 
stars. In this study, we will use a DF designed to resem- 
ble the old stellar disk. The exponential surface density 
has scale length R s — 0.33i?o> the exponential velocity 
dispersion is normalized to o-r(Rq) = 0.2vq and has scale 
length R(j = Rq. I also experimented with R a — 0.66i?o, 
corresponding to o\ oc E, and found very similar results. 

2.2.1. The Time Scales 

It turns out that the time t\, determining how smoothly 
the bar is switched on, is not very important for the out- 
come of the simulations. However, the total integration 
time under influence of the bar, i.e. w t-x — t%/2, has 
a significant impact on the details of the resulting ve- 
locity distributions. We will use ti — 2t\ and consider 
various values for t\ in units of the bar rotation period 
T b = 2vrAV 

2.2.2. The Solar Position 

The relative distance of the Sun from the bar is con- 
viently parametrized by the ratio Rolr/Rq-, where for 
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the power-law models 



Ro 



1 + 



1/(1-/3) 



(6) 



Another important parameter is the bar angle <f>, the 
angle by which the Sun is behind the bar's major axis 
(cf. Fig. @). 

2.2.3. The Rotation Curve 

The shape of the rotation curve, parametrized by the 
power (see equation [|a|), dictates the energy depen- 
dence of the orbital frequencies. We use a flat rotation 
curve (0 — 0) as default but will also consider slightly 
rising or falling rotation curves. 

Changing the local circular speed vq, or equivalently 
the frequency f2o, a t fixed R l.r/Ro corresponds to a 
change of the pattern speed fib by the same factor. Such 
a change is essentially a re-scaling of the model and does 
not require any additional simulations. 

2.24. The Bar Strength 

The strength of the bar is best measured in a way that is 
independent of such re-scaling. We will use as parameter 
the dimensionless quantity 



At 



Rb 

Rq 



(7) 



which is the ratio of the forces due to the bar's quadru- 
pole and the axisymmetric power-law background at 
galactocentric radius Ro on the bar's major axis. 

2.3. Scanning the Parameter Space 

Simulations have been performed for i? OLB /i?o between 
0.8 and 1.2 at steps of 0.05, nine values of at steps of 
0.05 between —0.2 and 0.2, and bar angles <p at 5° inter- 
vals between 10° and 50°. Additionally, the integration 
time and bar strength have been varied. In the figures 
below, unless stated otherwise, all the parameters but 
one are fixed to the default values listed in Table [3], in 
order to show the effect of one parameter alone. These 
default parameters correspond to an OLR less than 1 kpc 
inside the solar circle, a bar angle of 25° as discussed from 
studies of the gas motions, and a flat rotation curve. 

In order to ease a comparison with the observations, 
the velocity distributions arising from the simulations 
are displayed in the u and v, defined, respectively, as 
the velocity towards the Galactic center and in direction 
of rotation, both relative to the circular orbit passing 
through the position of the Sun in absence of any bar. 




0,0 



Fig. 2. — Simulated velocity distributions fo(u, v): variation with 
bar angle. The remaining parameters are fixed at their default 
values in Table hi 
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Fig. 3. — Simulated velocity distributions fo(u, v): variation with 
integration time t2 (T;, denotes the bar rotation period). The bar- 
growth time is set to ti = 0.5*2- 



Table 1 

Default values for the simulation parameters 



parameters to be varied 



default value 



shape of rotation curve 


13 





bar angle 




25° 


position w.r.t. OLR 




0.9 


bar strength 


Q 


0.01 


integration time 




4T b 


parameters kept fixed 


default value 


bar size 


Rb/Rcn 


0.8 


disk scale length 


Rs/Ro 


0.33 


local velocity dispersion 


a R (Ro)/v Q 


0.2 


or scale length 


Ra/Ro 


1 


bar growth time 


tl 


0.5 t 2 



2.3.1. Variations with Bar Angle 

Figure || plots fo(u, v) at intervals of 10° in cf> (note that 
the situation is bi-symmetric, so that cj> and (j> + 180° 
are identical) after the growth of a central bar. A bi- 
modality is clearly visible for bar angles in the first quad- 
rant. As judged from these simulations, the existence of 
the bi-modality in fo seems not to depend critically on 
the bar angle. 




Fig. 4. — Simulated velocity distributions fo(u, v): variation with 
Rol.r/ Ho- 

When increasing the bar angle, the depression be- 
tween the two modes, hereafter 'LSR mode' (at v ~ 0) 
and 'OLR mode' (at u, v < 0), is shifted towards higher 
v, while the OLR mode, which consists of stars origi- 
nating from inside the OLR, becomes more prominent. 
However, the latter effect is hard to compare quantita- 
tively with the observed fo(u,v), because the relative 
strength of the modes is likely to be affected both by the 
parameters of the DF and changes of the bar's pattern 
speed during the past ( Weinberg 1994 ). In u, the OLR 
mode ranges from about —0.6 to 0.2 in units of ^o. 

There is also a clear distortion of the LSR mode, 
which, for bar angles relevant to our position in the Milky 
Way, has the form of an extension to positive u velocities 
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-0.5 0.0 0.5 
u/v 



Fig. 5. — Simulated velocity distributions fo(u, v): variation with 
/3, parameterizing the shape of the rotation curve. 

at v comparable to that of the OLR mode. 

2.3.2. Variations with Integration Time 

Figure || plots fo(u,v) for three different choices of < 2 
(with ti = 0.5i2). Obviously, the longer integration 
times lead to more prominent resonant features, not only 
due to the OLR but also the resonance fit, — uJcf, + ldr, 
which is responsible for the ripple at vq ps 0.2. Clearly, 
the longer the action of non-axisymmetric forces, the 
stronger their imprints and the more detailed the result- 
ing structures in stellar velocities. Other effects of longer 
integration times are a slight shift of the peak v velocity 




-0.5 0.0 0.5 
u/v 



Fig. 6. — Simulated velocity distributions fo(u, v): variation with 
bar strength, parameterized by a. 

of the OLR mode towards lower values and a change in 
the structure of the LSR mode. However, the velocity 
of the minimum between the two modes hardly changes, 
justifying the use of t<i = ATi, as default. 

These simulations, assuming a constant pattern speed 
and a completely flat rotation curve without any sub- 
structures like spiral arms, will certainly create much 
more such fine resonant features in fo(u, v) than the more 
noisy and less constant force field of a real galaxy would 
produce. However, the locus of the resonance itself is not 
such a detail and we can safely assume that the v- velocity 
of the division line between the modes is not subject to 
artefacts due to the idealizations made. 

2.3.3. Variations with Radius 

Figure ^ plots fo(u,v) for five different choices of the 
observers distance from the OLR, where for R OLR /Ro = 
1.1 and 1.2, the integration time is taken to be t\ = 8T& 
in order to strengthen any possible resonant features. 
The bi-modality is clearly visible at 0.8 ;> Row./ Rq ;> 
1.05 - note, however, that this 'visibility' depends on 
the parameters of the DF and the bar angle <fi, cf. Fig. ||. 

For R OLR /R ^ 1.05, the OLR does not create a clear 
bi-modality (at least not at <j> = 25°) but a distortion in 
fo(u,v) in form of a ripple at roughly constant energy 
(curves of constant energy in the uv plane are circles 
centred on u — 0, v — —vq). For R l,h/Ro = 1-2, most 
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stars are on orbits between corotation and outer Lind- 
blad resonance. Despite this proximity to two funda- 
mental resonances, the resulting f (u,v) is remarkably 
regular: apart from the ripple caused by the OLR, it 
forms a nice elliptical distribution (bottom panel of Fig. 



1:1 sym 



1:1 asym 



2.3.4- Variations with Rotation- Curve Shape 

Figure || plots fo(u,v) for five different choices of /3, 
which determines the shape of the rotation curve (paj). 
Changing (3 affects the orbital frequencies, in the sense 
that they fall off more steeply with energy for a falling 
rotation curve (f3 < 0) than for a rising one (/? > 0) . For 
a fixed difference in frequencies, a steeper fall-off results 
in a smaller energy difference, which at fixed position 
yields a smaller velocity distance. Therefore, reducing [3 
squashes the structures of the velocity distribution in v. 
the extent in v of the OLR mode shrinks and its peak 
v velocity becomes less negative (the opposite happens 
when increasing (3) . Similarly, the distance in v between 
the OLR and higher-order resonances (responsible for 
the ripples at v > 0) decreases with decreasing (3 

2.3.5. Variations with Bar Strength 

Figure H plots fo(u,v) for three different bar strengths, 
parameterized by the force ratio a equation J?]]. Ap- 
parently, stronger bars lead to more pronounced OLR 
(as well as other resonant) features in fo(u,v): the OLR 
mode contains more stars and extends over a larger range 
in u velocity. 

The v velocity dividing the LSR and OLR modes is 
hardly changed. This is not what one would naively ex- 
pect from linear theory, in which the amplitude of the 
orbital changes is proportional to the amplitude of the 
perturbation. However, large-amplitude orbital changes 
are not suitably described by linear theory, but result, in 
the case of the OLR, in orbits that are detuned by the 
non-linearity and thus leave the near-resonant region of 
phase space. Thus, the OLR mode is made from orbits 
which are shifted by similar amplitudes largely indepen- 
dent of the strength of the bar perturbation. 

3. Relation to Closed Orbits 

When a bar is slowly grown, the initially circular or- 
bits will be mapped into (nearly) closed, but no longer 
circular, orbits in the barred potential. These are also 
the orbits on which gas is supposed to move, since en- 
counters are avoided. Thus, we expect most stars in a 
barred galaxy to move on nearly closed orbits. Conse- 
quently, the properties of the velocity distributions are 
better understood in light of the properties of the closed 
orbits that exist in the outer parts of a barred disk. 



ffi 1-5 



a; 



*r(2) 




0.5 



1 

L/(f R olr ) at 0-90" 



1.5 



Fig. 7. — Characteristic diagram (plot of R vs. L on the bar's 
minor axis) of the stable (solid) and unstable (dotted) closed orbits 
in the outer parts of the model with default settings. Orbits cor- 
responding to high-order resonances (n < —2) are not displayed. 
The dashed line corresponds to the circular orbits in the unbarred 
model with the square indicating the OLR. 

3.1. Closed Orbits in the Outer Disks of Barred 
Galaxy 

The only literature on this topic appears to be a review 
by Contopoulos & Grosb0l (1989), who consider closed 
orbits in the inner and outer parts of disks with bars or 
spiral structure. 

Outside of corotation, the main family of closed orbits 
are the near-circular xi orbits. At each integer value of 
the rotation number 



U>t/> 



(8) 



however, this family is modified by a resonance. At reso- 
nances with even n (as e.g. the OLR, which has n = —2), 
this modification takes the form of a gap with a distinct 
change in properties, such as shape, between the orbits 
on either side of the resonance. At resonances with odd 
n, two families of orbits, one symmetric and one anti- 
symmetric with respect to the bar's minor axis, branch 
off the main xi family. Immediately outside corotation, 
the denominator in equation (||) scans through all small 
negative numbers, resulting in a whole host of high-order 
resonances with — oo < n < — 2. Figure fj] shows the char- 
acteristics of the orbit families with n > —2 in the model 
with the default parameters of Section |2[ 

At n = —2 (OLR) we have the typical change, already 
mentioned in the introduction, from orbits inside OLR, 
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(a) x,(2) (b) xf(Z) 




Fig. 8. — Closed orbits of the families also displayed in Fig. 

The OLR is indicated by the bold circle. Solid and dotted 
curves refer to stable and unstable orbits, respectively. The xi(l) 
orbits extend to larger radii, while for the other four orbital families 
the innermost and outermost orbit is shown (note the difference 
in scale). The 1:1 resonant orbits (d & e) also exist as reflected 
versions, such that each family as a whole is bi-symmetric. 

which anti-align with the bar and are called xi(2) orbits 
by Contopoulos & Grosb0l, to orbits outside OLR, called 
xi(l) and aligning with the bar. The first subfamiliy, 
xi(2), extends to unstable orbits, called x*(2), at larger 
radii. Some orbits of these families are shown in the 
upper three panels of Figure ||. The innermost Xi(2) 
orbit clearly shows the effect of the n = —4 resonance. 

The resonance with n = —1, also called -1:1 or 
outer 1:1 resonance, creates a pair of orbit families (Fig. 
||d&e), each of which exists in two versions, which are 
reflection symmetric to each other. The family which is 
symmetric with respect to the bar's minor axis consists 
of stable orbits, while the antisymmetric family contains 
only unstable orbits. At large radii, both these family 
develop inner loops, which may penetrate right into the 
central bar. 



3.2. Closed Orbits and f 

All of the orbital subfamilies in Figure ^ pass through 
positions with R ~ R l.r, but there is no position (i?, 4>) 
which is visited by orbits from all five families. Thus, the 
accessibility by closed and nearly closed orbits changes 
between different positions in the outer Galaxy. Since 
each stable family (and non-closed orbits trapped around 
it) will create a distinct feature in the velocity distri- 
bution, we expect that fg(u,v) changes almost discon- 
tinuously between different positions, depending on the 
closed orbits families passing through that position. 

3.2.1. The Position Rolr/Rq = 0.9 and (f> = 25° 

Let us firts consider the case of the simulation with de- 
fault parameters, cf. the middle panels of Figs. ||and||. 
There are three closed orbits reaching this position: two 
xi(l) orbits and one x£(2) orbit. 

The first of the xi(l) orbits has (u,v) = (0.02, -0.015)u , 
while for the second (u,v) — (0.09, — 0.l)v . This latter 
orbit originates from very close to the OLR, resulting in 
a large pertubation, which enables its visit at Rq. Both 
these velocities lie within the LSR mode, but while the 
first is very close to the undisturbed circular orbit, the 
local velocity of the latter is related to the extension to- 
wards u > 0. Thus, we can understand these extensions 
of the LSR modes as caused by stars on nearly closed, 
hence strongly populated, and nearly OL-resonant, hence 
highly disturbed, orbits. 

The third orbit reaching (Rq , cf>) is unstable and has 
(u,v) — (—0.38, — 0.19)uo. This velocity is exactly in the 
valley between the two modes. Thus, it seems that this 
valley is caused by unstable, i.e. chaotic, orbits. 

3.2.2. Varying the Radius 

Let us now consider what happens if we change the dis- 
tance from the OLR. For 0.7 <J R lr/ Rq <J L, there is an 
unstable x*(2) orbit passing through, and we therefore 
expect a depression (a valley or a truncation) in fa. 

For R OIjR /Rq ;> 1 (the detailed value depends on (f>), 
there is a fundamental change, as then (i) no stable closed 
xi(l) orbits passes through, but (ii) a stable xi(2) in- 
stead of an unstable x£(2) orbit. These two changes are 
certainly related to the rather abrupt change of the mor- 
phology of fo(u,v) apparent in Figure ^, from a clear 
bi-modality to a mere density enhancement at roughly 
constant energy. 

3.2.3. Varying the Azimuth 

If we change the bar angle <fi at fixed R OLK /Ro, the fam- 
ilies of closed orbit passing through (Ro,(f>) change as 
well. In particular for R OLH /Ro ~ 1, the stable xi(2) 
orbits appear for <j) ;> 40°, which in Fig. ^ leads to a 
stronger OLR mode for ~ 90° than for 0^0°. 
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On the other hand, when increasing </>, the second 
stable xi(l) orbits disappear with the consequence that 
the extension of the LSR modes towards u > weakens 
as one moves from the bar's major axis to its minor axis 
in Fig. | 

3.2.4-. The Outer 1:1 Resonance 

Even though the stable 1:1 resonant closed orbits do not 
pass through positions near the OLR with \<j>\ <^ 40°, 
there is a clear resonant feature in the velocity distribu- 
tions at these positions. Hence, this feature, which takes 
the form of a density enhancement at certain velocities, 
must be due to related non-closed resonant orbits. The 
density enhancement is likely caused by the intersection 
of the various orbits of this family at 0° <J 4> <; 180°. 

4. QUANTIVYING THE OLR FEATURE 

In order to allow for a quantitative comparison of the 
OLR feature in the simulated fo(u, v) with the observa- 
tions in the next section, we need to use a well-defined 
quantity that is easy to measure and likely to be stable 
under deviations from the idealized conditions assumed 
in the simulations. Here, I will only use the w-velocity 
of the saddle point between the LSR and OLR mode. 
This velocity, hereafter denoted v lr, is essentially in- 
dependent of the details of the stellar DF and even the 
strength of the bar. 

4.1. An Axisymmetric Estimate 

We may first neglect the effect of the bar potential and 
estimate the velocities (at i?o) of the orbits that would 
become exactly resonant in the presence of a bar with 
pattern speed fib- For power-law potentials, the orbital 
frequencies are well approximated by the corresponding 
frequencies of the circular orbit with the same energy 
(Dehnen 1999a). Using this approximation and the rela- 
tions for power-law potentials (cf . Appendix B of Dehnen 
1999a) and neglecting terms 0(v 3 /vq), we may estimate 
that the local velocities of orbits that are in OLR satisfy 



u 2 _ _ 1 + /3 

— = w OLR = v - 

2v 1-/3 



1 - 



1 + V(l + /?)/2 



(9) 



Thus, the local velocities of resonant orbits form a 
parabola whose maximal v occurs at u = and de- 
pends on the distance from the OLR and the shape of 
the rotation curve. Note that the valleys between the 
two modes in the simulated fo(u,v) are actually nearly 
parabolic. They are, however, displaced from the u = 
axis and also in v. These displacements must be due to 
the quadrupole forces of the bar, neglected in the above 
estimate. 



Table 2 

Best-fit values for (a, b, c) in equation (ll 





a 


b 


c 


15° 


1.3549 


0.0761 


0.1362 


20° 


1.2686 


0.0642 


0.1120 


25° 


1.2003 


0.0526 


0.0892 


30° 


1.1424 


0.0406 


0.0711 


35° 


1.0895 


0.0298 


0.0538 


40° 


1.0420 


0.0200 


0.0423 


45° 


1.0012 


0.0103 


0.0316 


50° 


0.9653 


0.0012 


0.0238 



4.2. Quantifying the Local OLR Velocity 

The value for the velocity v OT , R is sensitive mainly to 
four parameters: the bar angle <f>, the relative distance 
Row./ ' Rq of the OLR, the normalization vq of the ro- 
tation curve, and its shape, parameterized by (3. The 
dependence on v is a simple scaling, while that on the 
other three parameters is less trivial. 

We might hope that the dependence on R OIjR /Ro and 
(3 is already largely described by u lr defined in equation 
(||) . Indeed, from plotting i; OL r at fixed bar angle versus 
5 lr, we find that to good accuracy the dependence on 
Rolr/ Ro and (3 is described by a linear function of the 
form 



Wolr ~ av OLH - (b + c[3)v . 



(10) 



For bar angles <j> G [15°, 50°], Table | lists the best-fit val- 
ues for a, b, and c obtained from fitting i> lr for R oljR /Ro 
in the range from 0.8 to 0.95. Note that the values for 
b and c are always small, i.e. v OLjR is largely given by 
o^olr- The RMS error made by this approximation is 
0.0035«Oj while the maximal error is 0.013i>o (occuring 
for iW/-Ro = 0.95, [3 = 0.2, <j) = 15°). 

5. Comparison with the Observed 
Velocity Distribution 

Figure ^ shows the distribution fo(u,v) for late-type 
stars in the solar neighborhood (note that the volume 
sampled by these stars corresponds to the size of the dot 
in Fig. [I] referring to a possible solar position). Here, 
the velocities u and v are relative to the local standard 
of rest (LSR) as measured by Dehnen & Binney (1998) 
from a sample of about 14000 stars in the Hipparcos cat- 
alogue (ESA 1997). This sample was constructed to be 
essentially magnitude limited in order to avoid any kine- 
matic biases. The distribution shown in Figure ^ has 
been inferred statistically ( Dehnen 1998| ) from the tan- 
gential velocities of about 6000 late-type main-sequence 
stars (B— V > 0.6 mag) and giants in Dehnen & Binney's 
sample. 
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Fig. 9. — The velocity distribution fo(u, v) inferred from Hippar- 
cos data for 3 527 main-sequence stars with B — V > 0.6 and 2 491 
mainly late-type non-main-scquenc e stars, high- velocity stars ex- 
cluded, in the solar neig hborhood jPchncn 199a ). The symbol 
indicates the solar velocity. Samples of early-type stars contribute 
almost exclusively to the low- velocity region (solid ellipse), which 
contains the most prominent moving groups. The region of inter- 
mediate velocities (broken ellipse) is mainly represented by late- 
type stars, of which ~ 15% fall into this region. Gray scales are 
linear in /o and the contours contain, from inside out, 2, 6, 12, 21, 
33, 50, 68, 80, 90, 95, 99, and 99.9 percent of all stars. The l-o 
uncertainty in the contour lines is between 0.3 and 3kms . 

The distribution fo(u,v) has been inferred by max- 
imizing the log-likelihood plus some penalty function 
which measures the roughness of / (see Dehnen 1998 
for more details). The latter ensures a smooth distribu- 
tion and suppresses the amplification of shot noise. Error 
estimates using the boot strap method (cf. Press et al. 
199^) indicate that the contours have an uncertainty of 
about 0.3-3 kms -1 . The two outermost contours are af- 
fected by the way the smoothing is introduced, and hence 
are less reliable. All features discussed in this section are 
of high significance. 

5.1. The OLR Induced Bi-Modality 

The distribution in Figure || shows a lot of structures. 
One of the most obvious is the bi-modality between 
a low-velocity component centred on the LSR motion 
(solid ellipse) and an intermediate-velocity component 
at v <^ —30 and predominantly negative u (broken el- 
lipse), which contains about a every sixth late-type star 
near the Sun. This clear bi-modality was, to the best 
of my knowledge, not known from pre-Hipparcos data; 
the corresponding mean outward motion of stars with 
v <; —30 was called 'it-anomaly'. The bi-modality is 
hardly present in samples of early-type stars, which al- 
most exclusively populate the moving groups (sub-peaks) 



in the low-velocity region ( 


Dehnen 1998 


; Chereul, Creze 


k Bienaymc 1998; Asiain 


et al. 1999 




Skuljan, Hearn- 


shaw & Cottrcll 1999). 









5.1.1. Is the Bi-Modality due to the OLR? 

The bi-modality present in the locally observed f (u 7 v) 
(Fig. ^|) is indeed very similar to those emerging from 
the simulations in Section ^|. Moreover, since we expect, 
from our previous knowledge of the bar and its estimated 
pattern speed, that -R lr ~ Rq (cf- Section Q), it is only 
natural to identify the observed bi-modality as the OLR 
feature of the Galactic bar. However, this may be pre- 
mature to do and we must first check if there are viable 
alternative explanations. 

First, can the intermediate-velocity mode at v w —45 
be the relic of a dispersed stellar (open) cluster or as- 
sociation (the standard explanation for the formation of 
moving groups)? The strongest argument against this 
hypothesis comes from the age of the participating stars, 
which from the absence of early-type stars with i nterme- 
diate velocities may be inferred to be > 8Gyr ( Dehnen 
1998 ). This high age, corresponding to ;> 40 orbital 
times, makes it very unlikely for any initial moving group 
to survive Galactic scattering processes. Moreover, the 
number of stars in the secondary mode is significantly 
higher than that in the low-velocity moving groups, such 
as the Hyades and Sirius streams. 

Second, could this secondary mode be the result of a 
merger with a globular cluster or satellite galaxy? This 
again can be ruled out, since such a scenario is highly un- 
likely to produce a feature with disk-like kinematics (a 
velocity deviating from the LSR motion by only ~ 15% 
instead of an expected 100%; a vertical velocity distribu- 
tion like that of low- velocity stars). Moreover, Raboud 
et al. (1998) have reported that the u-anomaly is caused 
by intermediate to high-metallicity stars, while a merger 
would involve rather metal poor stars. 

Third, could the bi-modality be caused by any other, 
possibly resonant, scattering process, e.g. due to spiral 
arm structure? While this possibility is harder to exclude 
than the first two, it seems rather unlikely. Spiral struc- 
ture should mainly affect stars with epicycles smaller 
than the inter-arm separation, while the bi-modality is 
created by stars with epicycles comparable or larger than 
3 kpc. The possibility of inner Lindblad resonant scatter- 
ing by some unknown agent corotating at ~ 25 kpc, e.g. 
a slowly rotating halo or a satellite galaxy (note, how- 
ever, that the magellanic clouds are too far out and have 
polar orbits, which disqualifies them as potential agents) 
was shown by Weinberg (1994) to be inconsistent with 
velocity dispersion data. 

So, we can conclude that outer Lindblad resonant 
scattering off the Galactic bar is presently the only viable 
explanation that can explain the bi-modality observed in 
fo(u,v) of late- type disk stars and its absence in early- 
type stars. 

Based on Fux's (1997) TV-body simulations of the 
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— v =180kms 
... v ^200 kms 
v„ = 220 km s 
— v„ = 240 kms 





30 
0[deg] 

Fig. 10. — The ratio Cl^/Slo that follows from the new constraint 
for given rotation curve parameters vo, f3, and bar angle <fr. The 
typical error in Q^/Qo due to the uncertainty in the measured dolr 
is about 1%, i.e. two small ticks on the y-axis. A small systematic 
error due to the v motion of the LSR is likely to lead to an over- 
estimation of Slb/f^O by up to a few % for <j> < 45° (see text). 



Milky Way, Raboud et al. (1998) and Fux (1999b) also 
attributed the it-anomaly to the influence of the Galac- 
tic baiQ, but were unable to investigate the parameter 
space and to achieve a resolution, both spatially and in 



vel ocity, comparable to that in the observed distribution 
(a general problem with TV-body s imu lations, see my re- 



ma ~ks at the beginning of Section O ) 



4 These authors, however, did not relate the effect to the OLR, 
but to the fact that the corresponding stars have Jacobi integrals 
Ej = E — Vl^L (the Hamiltonian in the corotating frame) that 
are just high enough to enable them to penetrate into the bar 
region itself. However, there are several hints pointing against 
this interpretation. First, the valley between the modes in the 
simulated f(v) is curved like lines of constant energy, while lines of 
constant E j are curved the other way around. Second, the change 
in the velocity seperating the modes when changing Aolr/Ao an d 
/3 cannot be described by a constant value of Ej corresponding 
to the value just allowing penetration into the bar. Moreover, the 
OLR does not occur at constant Ej (cf. Figs. 2, 3, and 4 of Fux, 
1999b), but roughly at constant E. 



5.1.2. Implications for the Inner Galaxy 

A quick comparison with the simulations shows that (i) 
the OLR must be inside the solar circle, i.e. .Rolr < Rq, 
and (ii) the bar angle <f> in the range 10° to 70°. We may 
even go further and make a quantitative comparison of 
the w-velocity 



Wolr = (-31 ± 3) kms 



,-i 



(11) 



of the valley between LSR and OLR mode with the cor- 
responding velocity measured in the simulated fo(u,v). 
In Figure flG, the ratio of the bar's pattern speed over the 
local circular frequency that is implied by the new con- 
straint via the approximation ([h]) is plotted for various 
possible choices for the local run of the rotation curve 
(as given by vq and (3) and the bar angle </>. 



The estimate in Fig. 10 is subject to some small sys- 
tematic errors. In particular, while in the simulations 
Wolr is defined relative to the circular speed vq in the 
absence of any bar, the value ([ll]) obtained from the ob- 
served fo(u,v) has been measured relative to the LSR, 
which shall deviate from an exactly circular orbit. From 
the simulations, it is simple to estimate the size of this 
deviation caused by the bar itself (for which it is small, 
especially if <j) ~ 45°), but this is not very meaningful as 
the LSR, being defined by low-eccentricity orbits, may 
well be affected by local spiral structure, not considered 
in the simulations. However, even if the azimuthal LSR 
motion as large as 10 kms -1 , it would introduce a sys- 
tematic error of only by 3% for Jlb/^Oj i-e. an amount 
that is still smaller than the uncertainty due to the un- 
known bar angle 4>- 

An IAU standard value of vq = 220 kms -1 and a flat 
rotation curve yield f2b between 1.8 and 1.9 times fio- 
When using the observed proper motion for the radio 
source Sgr A*, which is th ought to be associated with a 
supermassi ve black hole ( Eckart fc Gcnzel 1997 , Ghcz 



et al. 1998|) , to constrain the value of the lo cal cir cular 
frequ e ncy, one finds ilp ~ 28.5 kms -1 kpc -1 ( Ftcid et al. 
Backer fe Sramek 199£), and thus fib ~ 51 to 



1999 



54 km s 1 kpc 1 . 

See Dehnen (1999c) for a more dctaillcd analysis of 
the consequences for the value of f^b, when additionally 
to the proper motion of Sgr A* the terminal gas velocities 
are used to constrain the run of the rotation curve. The 
main result from that study is that the value for f2b de- 
pends only weakly on the assumed R$ and (f> (as long as 
they take reasonable values) and is thus rather narrowly 
constraint to be between 50 and 56 kms -1 kpc -1 . An- 
other result from that study is that corotation of the bar 
occurs at R cr /Rq « 0.5 to 0.6. This can be compared 
to the results from stellar- and hydro-dynamical model- 
ing of the inner Galaxy: Weiner & Sellwood (1999) and 
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Fux (1999a) report values in the same range, while En- 
glmaier & Gerhard (1999) obtain an upper limit of 0.5 
for R CR /R , i.e. slightly discrepant with the new value. 

5.2. More Similarities with the Simulations 

Apart from the low- vs. intermediate- velocity bi-modality, 
the observed velocity distribution in Figure || has also 
other features in common with the distributions obtained 
in the simulations of Section |2|. 

5.2.1. The Structure of the LSR Mode 

The structure of the main mode centred on the LSR, 
in particular the extension to u > at v rj — 0.1 «o ~ 
— 20kms~ 1 , is at least qualitatively similar, apart, of 
course, from the prominent moving groups. In the sim- 
ulations, this extension is caused by stars on orbits that 
are nearly closed and have large perturbative amplitudes. 
These orbits originate from near the OLR and therefore 
have large deviations from circularity such that they may 
nonetheless visit the solar neighborhood. 

However, since these orbits are nearly closed, they are 
liklely to be affected by spiral stucture and other local 
deviations from a smooth force field. As such effects have 
not been considered in the simulations of Section ^, it is 
not surprising that simulations and observations do not 
agree in detail. 

5.2.2. Features due to the outer 1 : 1 Resonance 

The ripples at (u, v) w (—80, —5) and (40, 5) are reminis- 
cent of the ripples, for instance in Fig. |§|, caused by the 
outer 1:1 resonance. This similarity persists even into 
details, which is more clearly recognizable when com- 
paring to simulations with longer integration times, for 
instance, the middle panel of Fig. 0. Here, the u < 
side of the ripple is at smaller v than the u > side, 
in agreement with the corresponding features in the ob- 
served fo(u, v). Moreover, there is a small excess of stars 
towards larger v at \u\ < OAvq, which may be compared 
to the 'bump' in the third and fourth contour (from out- 
side) at \u\ < 20kms _1 and v > 20kms _1 in Fig. ^. 
This bump is presumably caused by stars from outside 
the outer 1:1 resonance. 

6. Discussion and Conclusion 

The simulations presented in Section |^ showed how a ro- 
tating bar affects the distribution function of a surround- 
ing stellar disk. Since the non-axisymmetric component 
of the bar-induced forces falls off steeply with radius 
(cx r~ 4 for the quadrupole) reaching about 1% of the to- 
tal force at Ro, the influence of the bar is restricted to or- 
bits which are nearly in resonance with it. This influence 
is strongest for the outer Lindblad resonance (OLR), 



which occurs for orbits whose radial and azimuthal fre- 
quencies obey the relation f2b = w$ + hujR, where 17b is 
the bar's rotation rate (pattern speed). At orbits that are 
nearly in OLR, an otherwise smooth stellar distribution 
function becomes distorted. At any given position in the 
disk, this distortion is visible in the observable stellar ve- 
locity distribution f(u,v) if near-resonant orbits (i) pass 
through this position and (ii) are sufficiently populated. 
For a distribution function designed to resemble the old 
stellar disk of the Milky Way, the OLR leaves its imprint 
in the velocity distribution over a wide range of possible 
position in the disk. 

For positions (Ro, (f) orientated relative to the bar 
at angles 6 [10°, 70°] and situated outside the radius 
_/?olr where circular orbits arc in OLR, the resulting fea- 
ture in f(u,v) is a clear bi-modality: apart from the 
dominant stellar component with velocities similar to 
the local standard of rest (LSR), there is a secondary 
mode at u,« < 0. This OLR mode consists of stars 
with mean outward motion and slower rotation veloci- 
ties than the LSR. For R < R OLR , the OLR induced 
feature occurs at higher rotation velocities than the LSR 
and no longer takes the form of a clear bi-modality. As 
has been demonstrated in Section [| all these changes in 
the behaviour of fo(u,v) can be well understood from 
the properties of the closed orbits in the outer parts of a 
barred galaxy. 

Extensive simulations showed that the precise veloc- 
ity of the OLR-induced feature depends mainly on four 
parameters: the normalization and shape of the under- 
lying circular speed curve, v c (R), the observer's relative 
position with respect to the OLR, R OIjR /Rq, and their 
orientation relative to the bar, measured by the bar an- 
gle 4>. Note that the dependence on v c (R) is restricted 
to its run in the region visited by near-resonant slightly 
eccentric orbits, i.e. between about 0.7i?o and Rq for 
i?oLR = 0.9.Ro. These dependences are such that the 
(negative) v- velocity, v OLR , separating the OLR mode 
from the LSR mode of /o, increases (becomes smaller in 
modulus) with increasing bar angle </>, decreasing ratio 
Rolr/Rq, vo = Vc(Ro), and (3 = dlnw c /dlni?. 

The velocity distribution f(u, v) actually observed for 
late-type stars in the solar neighborhood (Fig. ||) shows 
indeed a secondary mode at u, v < 0, which is very simi- 
lar to the OLR modes emerging in the simulations. As I 
have argued in Section 5.1.1, there exists no satisfying ex- 
planation of this seeming anomaly other than being the 
feature induced by the OLR of the Galactic bar. Inter- 
estingly, the observed fo(u,v) has also other features in 
common with the simulated distributions. Firstly, exten- 
sions towards high u (inward moving) at slightly negative 
v may be related to a ridge in Fig. || at v ss — 20kms _1 
connecting the Hyades and Pleiades stream and reach- 
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ing at least to u w 60kms _1 . Secondly, ripples due to 
the outer 1:1 resonance seen in the simulated f(u,v) at 
v > and large \u\ may be related to the features in Fig. 
| at (u,v) = (-80,-5) and (40, 5) km s" 1 . 

First of all, this means that the structure of the local 
velocity distrubution for late-type stars provides a clear 
evidence for the existence of the Galactic bar, completely 
independent of any observation of the inner Galaxy itself. 
Moreover, this structure implies that the Sun is situated 
slightly outside of the radius -Rolr and at bar angles 
<j> between about 10° and 70°. One may even use the 
observed value, v OLR ., f° r the w-velocity seperating the 
OLR and LSR mode in order to deduce that the ratio 
between the bar's pattern speed Ob and the local circular 
speed fio is between 1.75 and 2. The uncertainty here is 
dominated by the uncertainty in the bar angle 4>, which 
according to other evidence must lie somewhere in the 
range 20 to 45 degrees. Combining this result with the 
value flo w 28.5kms _1 kpc -1 derived from the proper 
motion of the radio sour c e Sgr A* at the G alactic centre 
dBacker k Sramek 1999] |Reid et ol. 1999| ) yields Q h = 
(53 ±3)kms" 1 kpc" 1 . 

Clearly, the conclusion on the precise value of 
drawn on the basis of the simulations of Section || is 
subject to (presumably small) systematic errors, origi- 
nating from the simplifications made in the simulations 
and from the unknown v motion of the LSR. To improve 
on this, one needs to (i) use a more realistic model for the 
Galactic potential, (ii) compare not only v OLR of the sim- 
ulated fo(u,v) but also the feature of the 1:1 resonance 
(to reduce systematic errors due to the unknown LSR 
azimuthal motion), and (iii) combine this with hydro- 
dynamical simulations, which can be compared directly 
to the gas velocities observed in the entire inner Galaxy, 
constraining both the bar parameters and the rotation 
curve. Such simulations, in conjunction with the con- 
straint set by the proper motion of Sgr A*, are likely to 
determine the pattern speed, orientation, and structure 
of the bar as well as the rotation curve and hence mass 
distribution of the inner Galaxy. In such an analysis, the 
constraints set by the local fo(u, v) constitute important 
new ingredients as they largely reduce the freedom for 
the pattern speed f2t>, which is not very well constraint 
by modeling of the inner Galaxy as is apparent from the 
diverging results obtained in the past by different mod- 
ellers. 

Thus, the velocity distribution observed locally im- 
poses important new constraints for the structure of the 
inner Galaxy independent of and complementary to the 
photometry and kinematics, both stellar and gaseous, of 
that region itself. 
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